Zi spin liquids in S = 1/2 Heisenberg model on kagome lattice: projective symmetry 
group study of Schwinger-fermion mean-field states 
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Due to strong geometric frustration and quantum fluctuation, S = 1/2 quantum Heisenberg 
antiferromagnets on the kagome lattice has long been considered as an ideal platform to realize spin 
liquid (SL), a novel phase exhibiting fractionalized excitations without any symmetry breaking. 
A recent numerical study[l| of Heisenberg 5=1/2 kagome lattice model (HKLM) shows that in 
contrast to earlier results, the ground state is a singlet-gapped SL with signatures of Z2 topological 
order. Motivated by this numerical discovery, we use projective symmetry group to classify all 20 
possible Schwinger-fermion mean-field states of Z2 SLs on kagome lattice. Among them we found 
only one gapped Z2 SL (which we call ^[0, n]/3 state) in the neighborhood of [/(l)-Dirac SL state. 
Since its parent state, i.e. {/(l)-Dirac SL is foundQ] to be the lowest among many other candidate 
U(l) SLs including the uniform resonating- valence-bond states, we propose this Z2[0,ir]/3 state to 
be the numerically discovered SL ground state of HKLM. 

PACS numbers: 71.27.+a, 75.10.Kt 
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I. INTRODUCTION 

At zero temperature all degrees of freedom tend to 
freeze and usually a variety of different orders, such as 
superconductivity and magnetism, will develop in dif- 
ferent materials. However, in a quantum system with 
a large zero-point energy, one may expect a liquid-like 
ground state to exist even at T = 0. In a system 
consisting of localized quantum magnets, we call such 
a quantum-fluctuation-drivcn disordered ground state a 
quantum spin liquid (SL)2. It is an exotic phase with 
novel "fractionalized" excitations carrying only a frac- 
tion of the electron quantum number, e.g. spinons which 
carry spin but no charge. The internal structures of these 
SLs are so rich that they are beyond the description of 
Landau's symmetry breaking theory^ of conventional or- 
dered phases. Instead they are characterized by long- 
range quantum entanglement^ encoded in the ground 
state, which is coined "topological order"— & in contrast 
to the conventional symmetry-breaking order. 

Geometric frustration in a system of quantum magnets 
would lead to a huge degeneracy of classical ground state 
configurations. The quantum tunneling among these 
classical ground states provides a mechanism to realize 
quantum SLs. The quest for quantum SLs in frustrated 
magnets (for a recent review see Ref. @) has been pur- 
sued for decades. Among them the Heisenberg 5=1/2 
kagome lattice model (HKLM) 
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has long been thought as a promising candidate. Here 
< i,j > denotes i,j being a nearest neighbor pair. 
Experimental evidence of SL^— has been observed in 
ZnCu3(OH)gCl2 (called herbcrtsmithite) , a spin-half an- 
tifcrromagnet on the kagome lattice. Theoretically, in 



lack of an exact solution of the two-dimensional (2D) 
quantum Hamiltonian ([1} in the thermodynamic limit, 
in previous studies either a honeycomb valence bond 
crystal^— (HVBC) with an enlarged 6 x 6-site unit cell, 
or a gapless SlJ^ were proposed as the ground state of 
HKLM. However, recently an extensive density-matrix- 
renormalization-group (DMRG) studyi on HKLM re- 
veals the ground state of HKLM as a gapped SL, which 
substantially lowers the energy compared to HVBC. Be- 
sides, they also observe numerical signatures of Z2 topo- 
logical order in the SL state. 

Motivated by this important numerical discovery, we 
try to find out the nature of this gapped Z2 SL. Dif- 
ferent Z2 SLs on the kagome lattice have been previ- 
ously studied using Schwinger-boson representatio n 20 ' 21 . 
Here we propose the candidate states of symmetric Z2 
SLs on kagome lattice by Schwinger-fermion mean field 
approach 2 ^— . Following is the summary of our results. 
First we use projective symmetry groups (PSG) to clas- 
sify all 20 possible Schwinger-fermion mean-field ansatz 
of Z2 SLs which preserve all the symmetry of HKLM, 
as shown in TABLE |TJ We analyze these 20 states and 
rule out some obviously unfavorable states: e.g. gapless 
states, and those states whose 1st nearest neighbor (n.n.) 
mean-field amplitudes must vanish due to symmetry. 
Then we focus on those Z2 SLs in the neighborhood of 
the [/(l)-Dirac SB 2 .. In Ref. Hit is shown that £/(l)-Dirac 
SL has a significantly lower energy compared with other 
candidate U(l) SL states, such as the uniform resonating- 
valcnce-bond (RVB) state(or the U(l) SL-[0, 0] state in 
notation of Ref. 0). We find out that there is only one 
gapped Z2 SL, which we label as ^2[0,7r]/3, in the neigh- 
borhood of (or continuously connected to) Z7(l)-Dirac 
SL. Therefore we propose this Z2 [0, n]/3 state as a promis- 
ing candidate state for the ground state of HKLM. The 
mean-field ansatz of ^[0, tt]/3 state is shown in FIG-Hfb). 
Our work also provides guideline for choosing variational 
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FIG. 1: (color online) (a) kagome lattice and the elements 
of its symmetry group. Si,2 are the translation unit vec- 
tors, Ce denotes n/3 rotation around honeycomb center 
and cr represents mirror reflection along the dashed blue 
line. Here u a and up denote 1st and 2nd nearest neigh- 
bor (n.n.) mean-field bonds while u 1 and u~, represent two 
kinds of independent 3rd n.n. mean-field bonds, (b) Mean- 
field ansatz of ^[0, n]f3 state up to 2nd nearest neighbor. 
Colors in general denote the sign structure of mean-field 
bonds. Dashed lines denote 1st n.n. real hopping terms 

XiT, <l , ]>a (^jfLfj^ + h - c -)- red ones have V H = 1 and 
black ones have i/y — —1. Solid lines stand for 2nd n.n. 

hopping X2 E«y» a v ij(fLfj<* + h - c -) and singlet pairing 
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pVij(&2f} a f}p + h.c): again red ones have 
— 1. Here xi,2 and A2 
are real parameters after choosing a proper gauge. 



states in future numeric studies of SL ground state on 
kagome lattice. 



II. SCHWINGER-FERMION CONSTRUCTION 
OF SPIN LIQUIDS AND PROJECTIVE 
SYMMETRY GROUP (PSG) 

A. Schwinger-fermion construction of symmetric 
spin liquids 

In the Schwinger-fermion construction 2 ^—, we repre- 
sent a spin-1/2 operator at site i by fermionic spinons 
{fia, a=t,4}: 



(2) 



Heisenbcrg hamiltonian H 
represented as H = J2<ij> ~l J ij (fLfj<*f}pfjp + 
^flafiafjpfjp)- This construction enlarges the Hilbcrt 
space of the original spin system. To obtain the physical 
spin state from a mean- field state of /-spinons, we need 
to enforce the following onc-/-spinon-per-site constraint: 



1. 



fiafifiCaP — 0. 



(3) 



Mean- field parameters of symmetric SLs are Aije a p = 

-2(fiafj/3), Xijda/3 = mlafjf)), where e Q p is the com- 
pletely antisymmetric tensor. Both terms are invariant 



under global SU(2) spin rotations. After a Hubbard- 
Stratonovich transformation, the lagrangian of the spin 
system can be written as 

i <ij> 

- (417^ + h.c.)] + ]T a l Q {i)4r% (4) 



where two-component fermion notation tpi = (/if, fj, ) is 
introduced for reasons that will be explained shortly. We 
use r° to denote the 2x2 identity matrix and r 1,2 ' 3 are 



the three Pauli matrices, 
amplitudes: 

Ua = I 



Ua is a matrix of mean-field 



(5) 



a l (i) are the local lagrangian multipliers that enforce the 
constraints Eq.Q. 

In terms of ip 7 Schwinger-fermion representation has 
an explicit SU (2) gauge redundancy: a transformation 
ipi -> Wiipi, Uij -> WiUijW], Wi G SU(2) leaves the ac- 
tion invariant. This redundancy is originated from repre- 
sentation Eq.([2|): this local SU(2) transformation leaves 
the spin operators invariant and does not change physical 
Hilbert space. One can try to solve Eq.Q by mean-field 
(or saddle-point) approximation. At mean-field level, Uij 
and a are treated as complex numbers, and a l must be 
chosen such that constraints ((3|) are satisfied at the mean 
field level: (ipjT l ipi) = 0. The mean-field ansatz can be 
written as: 



H 



MF 



4a l Q T 1 ^. (6) 
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where we defined = | •/, , (',,. Under a local SU{2) 
gauge transformation —> Wi(i\j)Wj , but the physi- 
cal spin state described by the mean-field ansatz 
remains the same. By construction the mean-field ansatz 
does not break spin rotation symmetry, and the mean 
field solutions describe SL states if lattice symmetry is 
preserved. Different ansatz may be in different 

SL phases. The mathematical language to classify differ- 
ent SL phases is projective symmetry group (PSG)-. 



B. Projective symmetry group (PSG) classification 
of topological orders in spin liquids 

PSG characterizes the topological order in Schwinger- 
fermion representation: SLs described by different PSGs 
are different phases. It is defined as the collection of all 
combinations of symmetry group and SU (2) gauge trans- 
formations that leave mean-field ansatz invariant 
(as a l are determined self-consistently by these 
transformations also leave a l invariant). The invariance 
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of a mean- field ansatz under an element of PSG 

GjjU can be written as 

GuU({{i\j)}) = {{i\j)}, (7) 

u({(i\j)}) = {(i\j) = (u- l (i)\u- 1 (j))}, 

G u ({(i\j)}) = {(i\j) = G u (i)(i\j)G u (j)% 
G v (i) G SU(2). 

Here U G 5G is an element of symmetry group (SG) of 
the corresponding SL. In our case of symmetric SLs on 
the kagome lattice, we use (x, y, s) to label a site with 
sublatticc index s = u,v,w and 1,5 6 Z. Bravais unit 
vector are chosen as a\ = ax and 0,2 = §(i + V3y) as 
shown in FIG. [Ha). The symmetry group is generated by 
time reversal operation T, lattice translations along 
ai.2 directions, 7r/3 rotation C§ around honeycomb pla- 
quette center and the mirror reflection cr (for details sec 
Appendix |A")) . For example, if U = T\ is the translation 
along ai-direction in FigQJa), T\({x, y, s}) = {x+l,y, s}. 
Gjj is the gauge transformation associated with U such 
that GjjU leave invariant. Notice this condition 

([7]) allows us to generate all symmetry-related mean-field 
bonds from one by the following relation: 

(i\j) = G u (i){U- 1 (i)\U- 1 (j))Gl(j) (8) 

There is an important subgroup of PSG, the in- 
variant gauge group (IGG), which is composed of 
all the pure gauge transformations in PSG: IGG = 
{{Wi}\Wi(i\3)W} = (i\j),Wi G SU(2)}. In other words, 
Wi = G e (i) is the pure gauge transformation associated 
with identity element e G SG of the symmetry group. 
One can always choose a gauge in which the elements 
in IGG is site-independent. In this gauge, IGG can be 
the global Z2 transformations: {G e (i) = G e = ±r }, 
the global U(l) transformations: {G e (i) = e 10T ,6 G 
[0, 27t]}, or the global SU(2) transformations: {G e (i) = 
e i9fl ' ? ,e G (0, 27r],n G S 2 }, and wc term them as Z 2 , 
U(l) and SU(2) state respectively. 

The importance of IGG is that it controls the low en- 
ergy gauge fluctuations of the corresponding SL states. 
Beyond mean- field level, fluctuations of and a l need 
to be considered and the mean- field state may or may not 
be stable. The low energy effective theory is described 
by fermionic spinon band structure coupled with a dy- 
namical gauge field of IGG. For example, Z2 state with 
gapped spinon dispersion can be a stable phase because 
the low energy Z2 dynamical gauge field can be in the 
deconfincd phas o 29 ' 30 . 

Notice that the condition {G e (i) = G e = ±r } for 
a Z2 SL leads to a series of consistent conditions for the 
gauge transformations {Gjj{i)\U G SG}, as shown in Ap- 
pendix [A] Gauge inequivalent solutions of these condi- 
tions (|Al|) - (|ATT|) lead to different Z 2 SLs. Soon we will 
show that there are 20 Z2 SLs on the kagome lattice that 
can be realized by a Schwingcr-fcrmion mean-field ansatz 

am- 



III. Z 2 SPIN LIQUIDS ON THE KAGOME 
LATTICE AND Z 2 [0,7r]/3 STATE 

Following previous discussions, we use PSG to classify 
all possible 20 Z 2 SL states on kagome lattice in this 
section. As will be shown later, among them there 
is one gapped Z2 SL labeled as Z2[0,ir]f3 state in the 
neighborhood of £7(1)-Dirac SL. This Z 2 [0,tt]/3 SL state 
is the most promising candidate for the SL ground state 
of HKLM. 



A. PSG classification of Z2 spin liquids on kagome 
lattice 

Applying the condition G e (i) = G e = ±r° to kagome 
lattice with symmetry group described in Appendix 
\X[ we obtain a series of consistent conditions for the 
gauge transformation Gjj(i), i.e. conditions (|A4[) - (IA1I I) . 
Solving these conditions we classify all the 20 different 
Schwingcr-fermion mean-field states of Z2 SLs on kagome 
lattice, as summarized in TABLE [I] These 20 mean-field 
states correspond to different Z2 SL phases, which cannot 
be continuously tuned into each other without a phase 
transition. 

As discussed in Appendix IB 21 from PSG elements 
Gjj(i) one can obtain all other symmetry- related mean- 
field bonds from one using symmetry condition ([8]). 
Therefore we use u a = (0, 0, u|0, 0, u) to represent 
1st nearest neighbor (n.n.) mean-field bonds, up = 
(0, 1, w\0, 0, u) is the representative of 2nd n.n. mean- 
field bonds. There are two kinds of symmetry- unrelated 
3rd n.n. mean-field bonds, represented by u 7 = 
(1, 0, u|0, 0, u) and u 7 = (1, — 1, u\Q, 0, u). The symme- 
try conditions for these mean-field bonds are summarized 
in (|B13[) - (|B16[) . Besides, the on-site chemical potential 
terms A(i) (which guarantee the physical constraint ([3]) 
on the mean-field level) also satisfy symmetry conditions 
(|B12|) . We can show that A(x,y,s) = A s for these 20 
Z2 SL states. The symmetry-allowed mean-field ampli- 
tudes/bonds are also summarized in TABLE |TJ 

From TABLE[T]we can see there are 6 states, i.e. #7 — 
#12 that don't allow nonzero 1st n.n. mean-field ampli- 
tudes due to symmetry. Moreover, they cannot realize Z2 
SLs with up to 3rd n.n. mean-field amplitudes. Therefore 
they are unlikely to be the HKLM ground state. Ruling 
out these six Z2 SLs, we can see the other 14 Z 2 SL states 
fall into 4 classes. To be specific, they are continuously 
connected to different parent U(l) gaplcss SL states on 
kagome lattice. These parent U(l) SL states in general 
have the following mean-field ansatz 

Hu(i)sl = Xi v ^ULf 3a + h.c.) (9) 

<ij> 

where Vij = ±1 characterizes the sign structure of hop- 
ping terms with xi <= ^- Different parent U(l) SL 
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TABLE I: Mean-field ansatz of 20 possible Z2 SLs on a 
kagome lattice. In our notation of mean-field amplitudes 
{x, y, s\0, 0, u) = [x, y, s], this table summarizes all symmetry- 
allowed mean-field bonds up to 3rd n.n., i.e. 1st n.n. bond 
u a = [0,0, v], 2nd n.n. bond up = [0,1, w], 3rd n.n. bonds 
u-y = [1,0, u] and ii-y = [1,-1, u] as shown in FIG. QJa). A 3 
denote the on-site chemical potential terms which enforce the 
constraint (|11[) . r° is 2 x 2 identity matrix while r 1 ' 2,3 are 
three Pauli matrices, r 0,3 denote hopping while r 1 ' 2 denote 
pairing terms. means the corresponding mean-field ampli- 
tudes must vanish due to symmetry. Red color denotes the 
shortest mean-field bonds necessary to realize a Z 2 SL. In 
other words, the mean-field amplitudes with red color break 
the U(\) gauge redundancy down to Z2 through Higgs mech- 
anism. So in #3, #19 and #7 — #12 states a Z2 SL cannot be 
realized with up to 3rd n.n. mean-field amplitudes. Note that 
#15 state needs only 3rd n.n. bond u 7 to realize a Z2 SL (u 7 
not necessary) , while state #20 needs only u 7 to realize a Z2 
SL (it 7 not necessary) . Notice that when r\\2 = — 1 the mean- 
field ansatz (instead of the SL itself) will break translational 
symmetry and double the unit cell. There are six Z 2 SLs, 
i.e. #7 ' — #12 that don't allow any 1st n.n. mean-field bonds. 
Among the other 14 Z2 SLs with nonvanishing 1st n.n. mean- 
field bonds, only five Z 2 SL states, i.e. #1, #2, #5, #13, #15 
have gapped spinon spectra. #2 or z? 2 [0,7r]/3 state in neigh- 
borhood of [/(l)-Dirac SL is the most promising candidate of 
Z2 SL for the HKLM ground state. 



states are featured by the flux of /-spinon hopping phases 
around basic plaqucttc: honeycombs and triangles on the 
kagome lattice. 

The simplest example is the so-called uniform RVB 
state with = +1 for all 1st n.n. mean-field 
bonds. The hopping phase around any plaquette is 
1 = cxp[i0], and the corresponding flux is [0, 0] for [trian- 
gle, honeycomb] motifs. The 4 possible Z 2 spin liquids in 
the neighborhood^ of uniform RVB states (i.e. U(l) SL- 
[0, 0] state in Ref. 0) are classified in Appendix |DJ They 



WAYA (b)/oY(i>. 

X 7T ;< 7T X 7T > / 0X0 X > 
\7T/ VT/ 7 W / W w N / 

\!A!A!7 <,°.A .A_°/ 

FIG. 2: (color online) Mean-field ansatz of (a) 17(1) SL-[7r,7r] 
state and (b) U(l) SL-[7r,0] state, with 1st n.n. real hopping 

terms H M f = XiY, < i J>a ( u ij fLfj<* + h - c )- Colors again 
denote the sign structure of mean-field bonds: red dashed 
lines have Vi* = +1 and black dashed lines have Vij — — 1. 



arc #1, #5, #15, #13 in TABLE [Q and TABLEHH We la- 
bel them as Z 2 [0,0]A, Z 2 [0,0]B, Z 2 [0,0]C and Z 2 [0,0]D 
states. They all have gapped spectra of spinons. 

The ansatz of two other parent f (1) SLs are shown in 
FIG. [2] They both have 7r-flux piercing through a tri- 
angle basic plaquette. Following the above notations of 
hopping phase in [triangle, honeycomb] motifs, with ei- 
ther 7r-flux or 0-flux through the honeycomb plaquette, 
they are called U(l) SL-[7r,7r] state and U(l) SL-[7r, 0] 
state. There are three Z 2 SLs in the neighborhood of 
both U(l) SL states, i.e. #3, #17, #19 around U(l) SL- 
[tt, tt] state and #4, #18, #20 around £7(1) SL-[tt, 0] state. 
All these six Z 2 SLs have gaplcss spinon spectra, inher- 
ited from the two parent gapless U (1) SLs. To be precise, 
the spinon band structure of these six Z 2 SL states are 
featured by a doubly-degenerate flat band and a Dirac 
cone at Brillouin-zone center. This is in contrast to the 
numerically observed gap in two-spinon spectrum^, thus 
we can also rule out these 6 Z 2 SLs for the HKLM ground 
state. 

Another U(l) SL state is the so called C/(1)-Dirac SL 
or U(l) SL-[0, tt] state. Its mean-field ansatz is shown 
by the 1st n.n. bonds in FIG. [IJb). Clearly 7r-flux 
pierces through certain triangle plaquette with no flux 
through the honeycomb plaquette. According to varia- 
tional Monte Carlo studies^^, this [/(l)-Dirac SL have 
substantially lower energy compared to many other com- 
peting phases, including the uniform RVB state. There- 
fore we shall focus on those Z 2 SLs in the neighborhood 
of the [7(1)-Dirac SL in our search of the HKLM ground 
state. We need to mention that although unlikely, the 
four Z 2 SLs in the neighborhood of uniform RVB state, 
or U(l) SL-[0,0] state are potentially possible to be the 
HKLM ground state. 

In a previous study using PSG in Schwinger-boson 
representation^, it was shown that there are 8 differ- 
ent Schwinger-boson mean-field ansatz of Z 2 SLs on the 
kagome lattice which preserve all lattice symmetry. How- 
ever, these 8 Z 2 SLs may or may not preserve time- 
reversal symmetry. One can show that requiring all 
lattice symmetry and time-reversal symmetry, there are 
16 different Schwinger-boson Z 2 SLs on the kagome lat- 
tice. The relation between the 20 Z 2 SLs in Schwingcr- 
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fcrmion representation (see TABLE fl} and the 16 Z 2 SLs 
in Schwinger-boson representation are not clear. To clar- 
ify the relation between SL states in these two different 
representations, one can compare the neighboring (or- 
dered) phases of the SLs, e.g. by computing the vison 
quantum numbers^ of SL states. 



B. Z2[0,n]f3 state as a promising candidate for the 
HKLM ground state 

How to find those Z2 SLs in the neighborhood of (or 
continuously connected to) the [/(l)-Dirac SL? Naively, 
we expect the mean-field ansatz of these Z2 SLs can be 
obtained from that of [/(l)-Dirac SL by adding an in- 
finitesimal perturbation. To be specific, we require an 
infinitesimal spinon pairing term on top of the [/(l)-Dirac 
SL mean-field ansatz © or (|C1[) to break the IGG from 
U(l) to Z2 through Higgs mechanism. Mathematically, 
we need to find those Z 2 SL states whose PSG is a sub- 
group of the [/(l)-Dirac SL's PSG^i. Such Z 2 SL states 
are defined to be in the neighborhood of [/(l)-Dirac SL. 
Similar criterion applies to the neighboring Z 2 SL states 
of any parent U(l) or ST/ (2) SL state. 

We find out all four Z 2 SLs in the neighborhood 
of [/(l)-Dirac SLs in Appendix [Cj They are states 
#6, #2, #14, #16 in TABLE [fl labeled as Z 2 [0,7r]a, 
Z 2 [0,n]f3, Z 2 [0,7r]7 and Z 2 [0,7r](5 states respectively. 
Since the effective theory of a t/(l)-Dirac SL is an 8- 
componcnt Dirac fcrmion coupled with dynamical (7(1) 
gauge field£i22, we can find out all symmetry-allowed 
mass terms that can open up a gap in the Dirac-likc 
spinon spectrum. Following detailed calculations in Ap- 
pendix [C] we can see that among the four Z 2 SLs around 
the [/(l)-Dirac SL, only one state, i.e. Z 2 [0,ir]/3 (state 
#2 in TABLE [T] and [TT| can generate a mass gap in the 
spinon spectrum. In other 3 states the Dirac cone in 
spinon spectrum is protected by symmetry. The mean- 
field ansatz of Z 2 [0, tt)/3 SL state up to 2nd n.n. is shown 
in FIG.QIb): 

Hmf = E,(A 3 E a flfia + Ai/^/i + h.c.) (10) 

+Xl T, <ij>a "n >/;„/,„ + h.C.) + 
E«ij» "ij(X2 Ea fiafja + A 2 E Q/ 3 ^ P fUjfi + hX ') 

where e a ^ is the completely anti-symmetric tensor. We 
only list up to 2nd n.n. mean-field amplitudes because as 
shown in TABLE [J (see also Appendix [C]) , this Z 2 [0, 7r]/3 
state only needs 2nd n.n. pairing terms to realize a Z 2 SL. 
We can always choose a proper gauge so that mean-field 
parameters xi,2 and A 2 are all real. The sign structure 
of Vij = ±1 are shown in FIG. [ljb), with red denoting 
vij = +1 and other colors representing = — 1. As 
discussed in Appendix [C] the 2nd n.n. singlet-pairing 
term A 2 # not only break the U(l) gauge symmetry 
down to Z 2 , but also opens up a mass gap in the spinon 
spectrum. The on-site chemical potential Ai,3 are self- 



consitcntly determined by the following constraint: 

E,(Ea=t,;/l/-- 1 )= - (11) 

For further n.n. mean-field ansatz see discussions in Ap- 
pendix [Cj 

IV. CONCLUSION 

To summarize, motivated by the strong evidence of 
a Z 2 SL as the HKLM ground state in recent DMRG 
study-i, we classify all possible Z 2 SL states in Schwinger- 
fermion mean-field approach using PSG. We found 20 
different Schwingcr-fcrmion mean-field states of Z 2 SLs 
on kagome lattice, among which 6 states are unlikely due 
to vanishing 1st n.n. mean-field amplitude. In other 14 
Z 2 SLs only 5 possess a gapped spinon spectrum, which 
is observed in the DMRG result^. These five symmetric 
Z 2 SL states are all in the neighborhood of certain parent 
U(l) gapless SLs. To be precise, four are in the neigh- 
borhood of gapless uniform RVB (or U(l) SL-[0, 0]) state, 
while the other one, i.e. Z 2 [0, ir](3 is in the neighborhood 
of gapless [/(l)-Dirac SL (or U(l) SL-[0,tt]) state. Pre- 
vious variational Monte Carlo study^ showed that gap- 
less Z7(l)-Dirac SL has a substantially lower energy in 
comparison to the uniform RVB state. This suggests Z 2 
SLs in the neighborhood of Z7(l)-Dirac SL should have 
lower energy compared to those in the neighborhood of 
uniform RVB state. Therefore we propose this Z 2 [0,ir](3 
state with mean- field ansatz (|1Q[) shown in FIG. QJb) as 
the HKLM ground state numerically detected in Ref. [l|. 
Our work provides important insight for future numeric 
study, e.g. variational Monte Carlo study of Gutzwiller 
projected wavefunctions. 
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Appendix A: Symmetry group of kagome lattice and 
algebra conditions for Z2 spin liquids 

As shown in FIG.QJa), we label the three lattice sites in 
each unit cell with sublattice index {s = u, v, w}. Choos- 
ing Bravais unit vector as a\ = ax and a 2 = ^(x + v3j/), 
the positions of the three atoms in a unit cell labeled by 
indices i = (x, y, s) are 

f(x,y,u) = (x+ 5)01 + (y + 5)02, (Al) 
r(x, y, v) = (x + |)oi + ya 2 , 
f(x, y, w) = xdi + (y + i)a 2 . 
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The symmetry group of such a two-dimensional kagome 
lattice is generated by the following operations 

Ti: (x,y,s)^(x + l,y,s); (A2) 
T 2 : (x,y,a) -> (x,y + l,s); 
a : (x,y,u) (y,x,u), 
(x,y,v) (y,x,w), 
(x,y,w) (y,x,v); 
G 6 : (x,y,u)^ (-y-l,x + y + l,v), 
(x,y,v) (-y,x + y,w), 
(x, y, w) (-y - l,a; + y, u). 

together with time reversal T. 

The symmetry group of a kagome lattice is defined by 
the following algebraic relations between its generators: 

T 2 = a 2 = (G 6 ) 6 = e, (A3) 
g- 1 T- 1 gT=e, V g = T 1 , 2) o-,C 6 , 





^f^aTi = 


e. 




l T^rrT 2 = 


e. 


cr 


l T 2 l aT x = 


e. 




1 T 2 ^ 1 C 6 Ti -- 


= e i 




T^TiCeTi 


= e 


cr 




e. 



where e stands for the identity element in the symmetry 
group. Therefore the consistent conditions for a generic 
Z 2 PSGs on a kagome lattice is written as 



[G T (i)} 2 = Vtt°, 




(A4) 


G,r(tr(i))G a (i) = 


-V.r°, 




(A5) 


G^ i ( i )G^( i )G Tl (i)G T (T 1 - 


" 1 «) = 


= r] Tl TT°, 


(A6) 




" 1 «) = 


-- r]T 2 TT°, 


(A7) 


gUi)gU*)gAi)g t (<t- 


■ 1 (0) = 


VctTT , 


(A8) 


gI^)gU^G Cc ^)Gt(C 6 


_1 (0) = 




(A9) 



G^ 2 (T 1 - 1 (z))G^(z)Gt 2 ( 1 )G Ti (T 2 - 1 ( 1 )) = ?7 12 r (A10) 

G C6 (G 6 - 1 (i))G Ce (G 6 - 2 (i))G Ca (G 6 3 ( l ))G C8 (G 6 2 (i)) • 
Gc 6 (G 6 2 ( l ))Gc 6 (G 6 ( l ))G C6 ( i ) = r?c 6 r°, (All) 
Gi(T 2 - 1 ( l ))G^ 2 ( l )G (T ( l )G Tl ( C r( l )) = VTi t , (A12) 

Gt(Tf 1 ( l ))G^ i ( l )G CT ( l )G T2 (cr( l )) = ^r , (A13) 
Gt(G 6 (i))Gc 6 (G 6 (i))G CT ( l )G C8 (cr( l )) = ^r^AW) 
GL 6 (r 2 - 1 ( l ))G^ 2 (*)Gc 6 (*)G Tl (G 6 - 1 (z)) = r/ C6Tl ^15) 
G^(T 2 - 1 r 1 ( l ))G; 2 (T 1 (z))G Tl (TxW) • 

G^Gt^" 1 ^)) = r,c 6T2 r . (A16) 

for any lattice site i = (x, y, s). Here all 77s are Z 2 integers 
characterizing different SLs: different (gauge inequiva- 
lent) choices of these Z 2 integers (different Z 2 PSGs) cor- 
respond to different Z 2 SLs. Notice that under a local 



guagc transformation W(i) £ SU(2) the PSG element 
Gjj{i) transforms as 

GuW^W^GuiiWHu- 1 ^)) (A17) 

Appendix B: Classification of all Z2 spin liquids on 
kagome lattice 

1. Classification of Z2 algebraic PSGs on kagome 
lattice 

In this section we classify all possible Z 2 spin liquids 
on a kagome lattice. Mathematically we need to find out 
all gauge-inequivalent solutions of algebraic conditions 
(|A4ll - (|Al5|) for Z 2 PSGs. 

First from condition (|A10[) we can always choose a 
proper gauge so that 

G Tl (x,y,s)=r, y 12 T°, G T2 (x,y,s) = r°. (Bl) 

From (|A12|) and (| A13|) we can see Ga-(x,y,s) = 
7 7ctTi 7 7<tT 2 7 7i2$<t(s)- Condition (|A5[) further determine 
rjcrTj = VrTT 2 an d therefore we have 

Gcr(x,y,s) = rfj^r^g^s) 

where SU{2) matrices ga-(s) satisfy 

9*{w)g a {v) = [ g<T {u)] 2 = r 1(r r Q (B2) 

Notice that we can always choose a proper global 
Z 2 gauge on Gt x (x, y, s) (which doesn't change the 
mean- field ansatz) so that rjc 6 T 2 — 1 m (| A16|) . 
From (|A15|) and (|A16|) it's straightforward to show 
that G Ce (x,y,u/v) = Vo^^ii^ 2 9c e (u/v) and 

G ( \ x-\-y x-\-y-\-xy-\-x(x+l) /2 / \ TJ _, 

c e (x, y,w) = V c „T,Vi2 9c e {w). It's con- 

dition (|A14[) that determines ?7c 6 Ti = i}trT-J]vi an( i finally 
we have 

Gc 6 {x,y,u/v) =rf crT v i r]l 2 2 g Cs {u/v), 
G Cli {x,y,w) = (ri 1 2<qvT 1 Y +v ri 12 2 gc 6 (w). 
where SU(2) matrices gc 6 ( s ) satisfy 

[gc e {w)g Ce (v)gc e (u)] 2 = Vl2VC e T° , (B3) 

[.9<t(w)5c 6 (w)] 2 = gAw)gc e {v)g,r{u)gc e {u) = r\„r\ aCf ^. 

(B4) 

according to (|A11|) and (|A14|) . 

Now through a gauge transformation W(x, y, s) = 
rfrT we can ^ X VcrTi,2 — 1 an d the PSG elements be- 
come 

G<r(x, y, s) = rfZ%g<r(s); (B5) 

x x(x + l) 

Gc 6 (x,y,u/v) = rf x l 2 g Ce (u/v) , 
G Ce (x,y,w) =nll +x+v+ 2 gc s (w)- (B6) 
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II 

# 


?7l2 


g<r(u) 




3<t(™) 




5c 6 (» 


gc 6 (w) 


Label 


1 


+ 1 


u 

T 


u 

T 


u 

r 


u 

T 


u 

r 


u 

r 


Z 2 [0, OJA 


2 


— 1 


Ji 
T 


r 


_U 
r 


_U 

r 


_U 
T 


_u 
r 


Z 2 [0,7rJ£f 


6 


+ 1 


_u 

T 


T 


„u 
r 


„U 

T 


— r 


ir 




4 


— 1 


u 

T 


u 

T 


u 

T 


U 

T 


u 

— T 


• l 

IT 


v r nl A 

Z 2 [ty, 0\A 


r 




+ 1 


Ji 
T 


_U 
T 


_U 
r 


ir 


ir 


ir 


v rn ni n 
Z2 [0, \J\B 


6 


— 1 


U 

T 


U 

T 


u 

r 


j a 
ir 


■ a 

IT 


• a 

IT 


Z2 [0, 7r]a 


/ 


+ 1 


IT 


,,U 

r 


— T 


,,U 
r 


ir 


_u 

r 




o 



— 1 


lr 


T 


— T 


_U 

T 


lr 


_U 
r 




9 


+ 1 


• i 

IT 


u 

r 


u 

— r 


u 

r 


• i 

— IT 


• i 

ir 




1 n 
10 


— 1 


it 


T 


— T 


,,U 

T 


— IT 


ir 




11 


1 1 
+ 1 


1T 


_u 
r 


JJ 

— T 


ir 


— IT 


ir 




12 


— 1 


■ i 

it 


u 

r 


(J 

— T 


■ 3 

IT 


■ 2 

— IT 


• a 
ir 




to 


1 1 
+ 1 


it 


lr 


ir 


ir 


IT 


ir 


Z 2 [U, 0JL> 


1 /I 

14 


— 1 


it 


it 


ir 


ir 


IT 


ir 


Z 2 [0,7Tj7 


i r 

to 


I i 
+ 1 


IT 


IT 


ir 


„U 

T 


,u 
r 


_u 

T 


Z 2 [0, OJC 


ID 


1 

— 1 


1 T 


IT 


1 T 


_U 
T 


_U 

T 


_(J 
T 


Z 2 [U, 7TJO 


17 


+ 1 


ir a 


ir y 


ir 3 


r" 


r ,J 


ir 1 


^ 2 [7T,7r]B 


18 


-1 


ir ;! 


ir 3 


ir 3 


r° 


r u 


ir 1 


Z 2 [tt,0]B 


19 


+ 1 


ir a 


ir 3 


ir 3 


ir 3 


-ir 3 


it' 


Z 2 [n,n]C 


20 


-1 


ir ;j 


ir 3 


ir 3 




-ir ;l 


it 2 


Z2VAC 



TABLE II: A summary of all 20 gauge-inequivalent PSG's 
with Gt(x,i/, s) = ir 1 on the kagome lattice. Notice that 
there is a free Z2 integer 7712 = ±1 in other PSG elements 
(|B"TT) , (fB5l) and (fB6l) . They correspond to 20 different Z 2 spin 
liquids on the kagome lattice. 



According to (|A4[) . (|A6|) and (|A7|) we can see that 
Gt(x,u,s) = rj^ iT i^ T gT{s). (|A9|) and ()A8j) further de- 
termines Ty^T = t?t 2 t = 1 and by choosing a proper 
gauge we have 



G> 



r 

^x,y,s) = g T (s) = < 

IT 



= 1. 
?]T = -1. 



(B7) 



which satisfy 



gcr{u)g T (u) 
gv{v)g T {w) 
g<r{w)g T (v) 
9c 6 {u)g T {w) 
9cA v )9t{u) ■■ 
9c 6 (w)9t(v) ■- 



-- ^<tT5t(u).9 ct (w), 

= ?7 ct t.9t(w)5<t(«) i 
: r] rT T9T(w)g rT (w); 

Vc e T9T{u)gc 6 (u), 
■ Vc 6 TgT(v)gc 6 {v) 7 
Vc 6 TgT(w)g Ce {w). 



(B8) 



(B9) 



according to (|A9[) and 

In the following we find out all the gauge-inequivalent 
solutions of SU(2) matrices gT,a-,c e ( s ) satisfying the 
above conditions. They are summarized in TABLE . 

(I) g T (s) = t° and therefore i] T = r\ aT = Vc e T = 1: 
Conditions (|B8j) and (|B9j) are automatically satisfied, 
(i) V, = 1: 

Notice that under a global gauge transformation 
W(x,y, s) = W s G SU{2) the PSG elements transform 



as 

3 CT (w) -> W w g„{w)W$\ 
gc 6 {u) -> W„5c 6 (w)W+, 
3c 6 (w) -> W^5c 6 (u)W^, 

Thus from (|B2p and (|B4|) we can always have ga-(s) = t° 
and gc B (u) = r°, gc 6 («) = r/o-c 6 T by choosing a proper 
gauge. 

(A) 77 ctC6 = ?7i2??c 6 = 1: 

from (|B3[) we have gc 6 ) = T ° ■ 

(B) Tla-Ce = Vl2VC 6 = ~1: 

from (|B3|) we have gc 6 (w) = ir 3 by gauge fixing, 
(ii) = -1: 

from (|B2p we have ga-(v) = —g a {w) = t° and gcr(u) = 
ir 3 by gauge fixing. Also from (|B4[) we can choose a 
gauge so that gc 6 {u) = t° and gc 6 {v) = -^Vo-c 6 t 3 . 

(A) i^ce = -1: 

In this case (|B4[) requires gc 6 ( w ) — T ° an d thus 
Vi2Vc e = -1 according to (|B3|) . 

(B) ??(TC6 - 1: 

(a) r) 12 ric 6 = -1: 

Now from (|B4[) and (|B3|) we have gc 6 { w ) = 1 7 " 1 by 
gauge fixing. 

(b) 7? 12 ?/c 6 = 1: 

by (|B4[) and (|B3| we must have gc 6 ( w ) = 11-3 ■ 
To summarize there are 2 x (2 + 3) = 10 different 
algebraic PSGs with t]t = 1 and <?t(s) = r°. 

(II) <?t( s ) = 17-1 and ryx = — 1: 

(i) V, = 1: 

According to jB2|and (|B8|) . by choosing a proper 
gauge we can have <7 CT (s) = r° and r\ a T — 1. From 

dB3j) and (|B4]l wc also have [gceM] = 9c 6 {v)gc 6 ( u ) = 

V<tC 6 T° = Vl2VC 6 T°- 

(A) m2Vc 6 = V<rC e = 1: 

From (|B9[) , (|B3|) and (|B4|) , by choosing gauge we have 
gc 6 (s) = r° and t] C ' 6 t = 1. 

(B) ??i2?/c 6 = ??<tC 6 = -1: 

(a) ?7c 6 t = 1: 

In this case we have gc 6 ( u ) = ~9c 6 ( v ) — T ° and 
gc 6 { w ) = i 1 " 1 by choosing a proper gauge. 

(b) ??c 6 t = -1: 

In this case we can have gc 6 i s ) = 11-3 by choosing a 
proper gauge. 

(ii) Tla- = -1: 
(A) i^t = 1: 

From (|B8[) and (|B2p we have g CT (u) = ir 1 and gcr(v) = 
—gcr(w) = t° by proper gauge fixing. Also from (|B4|) 

we know [gc 6 {w)] 2 = -i^c 6 T ° and gc 6 {u)gc 6 {v) = 

-i^o-CsT 1 . 

(a) T)<rC 6 = -1: 
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from jB9| and jB4|, (JB3]) it's clear that t] C : 6 t = L 
9c B (u) = 9c 6 ( w ) = T ° and 9c e {v) = ir 1 through gauge 
fixing. Also we have f?i2f?c 6 — — 1- 

(b) 77 ctC6 = 1: 

(bl) r]c 6 T = 1: 

In this case r\n"r\Cf, = 1, and we can always choose a 
proper gauge so that gc 6 (u) = t°, gc 6 ( w ) = -9c e (v) = 
it 1 . 

(b2) VCeT = -1: 

In this case ?/i2??c 6 — ~L an d we can always choose a 
proper gauge so that gc e {v) = - it 2 , gc 6 ( u ) = 5c 6 (w) = 
ir 3 . 

(B) ri oT = -\: 

Conditions (|B8|) and (|B2|) assert that g„(s) = ir 3 by 
proper gauge choosing. 

(a) VcrCe = -1: 

In this case from (|B4j) we know gc e ( w ) — i 1 " 3 : hence 
Vc 6 t = — 1. Then we can always choose a gauge so that 
gc 6 (u) = gc 6 (v) = ir 3 and so r) 12 r) Ce = -1 from ([B"3]) . 

(b) r; CT c 6 = 1: 
(bl) r/c 6 r = 1: 

In this case from (|B8D . (IB4I) we have gc e ( u ) — 9c s (v) — 
r° by a proper gauge choice. Meanwhile conditions 
(|B3|) and (|B4[) become [gc 6 ( u ')] = Vi2VC 6 T° an d 
[ir 3 5 c 8 H] 2 --r°. 

(b.1.1) Time's = 1: 

here we have gc 6 ( w ) = r °- 

(b.1.2) TJi^Ca = -1: 

here we have gc e ( w ) = it • 

(b2) r /Cer = -1: 

In this case from (|B8|) and (|B4[) we can always choose 
a proper gauge so that gc 6 ( u ) = ~gc e ( v ) = i 1 " 3 - We also 
have gc 6 ( w ) = xt2 an d Vi2Vc e — — 1 from (|B3j) . 

To summarize there are 2 x (3 + 7) =20 different 
algebraic PSGs with rjx = — 1 and gr(s) = it 1 . 

So in summary we have 10 + 20 = 30 different Z-i 
algebraic PSGs satisfying conditions (|A4j) - (|A16|) . Among 
them there are at most 20 solutions that can be realized 
by a mean-field ansatz, since those PSGs with 5t(s) = 
t° would require all mean-field bonds to vanish due to 
(jBlip . As a result there are 20 different Z 2 spin liquids 
on a kagome lattice. 



2. Symmetry conditions on mean-field anstaz 

Let's denote the mean-field bonds connecting sites 
(0,0, u) and (x,y,s) as [x, y, s] = (x, y,s\0, 0,u). Using 
([8|) we can generate any other mean-field bonds through 
symmetry operations (such as translations Gt x 2 2~i,2 and 
mirror reflection G CT cr) from [x, y, s]. However these 
mean-field bonds cannot be chosen arbitrarily since they 
possess symmetry relation ([5]): 

= Gu(i) {U-Hm^U)) G f uU) (B10) 



where U is any element in the symmetry group. Notice 
that for time reversal T we have 

G T m\j)GUj) = -m (Bn) 

We summarize these symmetry conditions on the mean- 
field bonds here: 

(i) For s = u 

T: g T [x,y,u]g} r = -[x,y,u], 
Tf T^ x ct : [x, — x, u] — > [x, — x, u]', 

cr : [x, x, u] — > [x, x, u]. 

(ii) For s = v 

T: g T [x,y,v]gl, = [x,y,v], 

T*-*>T»- 1 *C 6 - 1 : [1 - 2y,y,v] -+[1- 2y,y,v] f - 

(iii) For s = w 

T: g T [x,y,w}gl r = [x,y,w\, 
T z-i T 2-2x a C 6 : [x,l- 2x,w] -)■ [x,l - 2x,w}\ 
T* +1 (tCq- 2 : [x,0,w] -> [x,0,w]l. 

Now let's consider several simplest examples. At first, 
on-site chemical potential terms A(x,y,s) = A s satisfy 
the following consistent conditions: 



T 1 hgT 1 = -A, 




g„(u)k u gl(u) = 


A u , 


g rT (v)A w gl(v) = 


A„, 


g (T {w)A v g} T {w) = 


A w ; 


gc 6 (u)A w gl B (u) = 


= A u , 


gc 6 (v)A u gl 6 {v) = 


-- A v , 


gc 6 (w)A v gl e (w) -- 


= A w . 



In fact in all 20 Z2 spin on a kagome lattice we all have 
A u = A v = A w = A s with a proper gauge choice. 

All the 1st n.n. mean-field bonds can be generated 
from u a = [0, 0, v] . For a generic Z 2 spin liquid with PSG 
elements G T (x,y,s) = it 1 and (|BT) (|B5)l (|B6|l . the bond 
u a = [0,0, v] satisfies the following consistent conditions: 

r 1 M ce r 1 = —u a , (B13) 
go-(u)gc 6 (u)gc 6 (w)u a gl B (v)gl 6 (w)gi(v) = u+ . 

It follows immediately that for six Z2 spin liquids, 
i.e. #7 - 12 in TABLE El all n.n. m ean-fi eld bonds must 
vanish since u a = as required by (|B13[) . Therefore it's 
unlikely that the Z2 spin liquid realized in kagome Hub- 
bard model would be one of these 6 states. In the follow- 
ing we study the rest 14 Z2 spin liquids on the kagome 
lattice. 
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All 2nd n.n. mean-field bonds can be generated from the above mean-field ansatz (|C4[) is 
up = [0,1, w] which satisfies the following symmetry con- 
ditions 



-up, 



(B14) 



9v{u)gcA u ) u f39c^ v )9i{w) = Ma- 



There are two kinds of 3rd n.n. mean-field bonds: the 
first kind can all be generated by it 7 = [1,0, u] which 
satisfies 

t^-u^t 1 = — « 7 , (B15) 
9c e (u)gc e (v)gc 6 (w)u^[gc 6 (u)gc 6 (v)gc e (w)] = u\. 

the second kind can all be generated by u 7 = [1, — l,u] 
which satisfies 



ff 1T (u){t 7 5i(u) = uL 



9C 6 {u)gc 6 (w)gc 6 (v)u>y [gc 6 (u)gc 6 (w)gc 6 (v)] ^ = T)12U 



G T (x, y, v) = G T {x, y, w) = -G T (x, y, u) = gr, 
grT 3 g f T = -r 3 ; 
Gt 2 (x, y, s) = g T2 , g T2 T 3 g^ 2 = t 3 ; 
G Tl (x, y, v) = G Tl (x, y, w) = -G Tl (x, y, u) 

= (-ir+y gTl , g T ,T 3 g Tl =r 3 ; 

G„{x,y,v)=G <T {x,y,w) = (-l)*+»+ 1 G tr (x,y,u) 

= { _ 1) ( X +y)(x + v + l)/2 g<Tj 5(tT 3^ = r 3. 

, \ / , g(g + l)-H/(«-l) 

Gc B {x,y,u) = (-1) ^ 5Ce; 
G Cs {x,y,v) = -{-!) ^ 5C( , ; 



Gc 6 ( x ^y^ w ) = 



j;(x-l) + B(i;-l) 



9Ce T i 



9c 6 T 3 gi CB = 



(C5) 



(B16) so ^ na ^ ^ ne m ean-field ansatz satisfy (|SJ). 



2. Classification of Z2 spin liquids around (7(1) 
SL-[0,tt] state 



Appendix C: spin liquids in the neighborhood of 

17(1) SL-[0,tt] state 

1. Mean-field ansatz of (7(1) SL-[0,vr] state 

Following ST/ (2) Schwinger fermion formulation with 
= {fi-\,fl^) T , we focus on those Z 2 spin liquids (SLs) 
in the neighborhood of U(l) SL-[0,7r] state with the fol- 
lowing mean-field ansatz: 

(x, y, u[x, y, v) = -(x, y, u[x, y, w) = (~1) x xt 3 , (CI 
(x + l,y,w[x,y,u) = (x,y + l,v\x,y,u) = -(x,y,v\x,y,w) 
= (x + 1, y - 1, w\x, y, v) = \t 3 . 

where \ is a rea l hopping parameter. We define mean- 
field bonds (x, y, s\x' , y' , s') in the following way 



Plugging (|C5|) into algebraic consistent conditions 
V4|) - (|A15[) yields four algebraic solutions of Z 2 PSGs 
around the U (1) SL-[0, 7r] state. Choosing a proper gauge 
they all satisfy 



9t = m, gr x = gr 2 = r , 
ryr = r)\i = flc^ = — 1, 

VT 1:2 T = TjcrT = T]c 6 T 2 = 1- 



(C6) 



The four Z2 PSGs near the Z7(l) SL-[0, 7r] state are fea- 
tured by 



9a- = 9Cs = T< \ 



(#6) Z 2 [0,7r]a 

9<r = Va-T = 1, 

VtC 6 = VC e T = -?7C 6 = ~1; 



(C2) 



(#2) Z 2 [0,n](3: g, 

9<r = Va-T = 1, 
?7<rC 6 = 7 ?C 6 T = -VG 



T°, 9C e = IT 3 , 



(C7) 



(C8) 



For convenience of later calculation we implement the 
following gauge transformation 

(C3) 

and the original mean-field ansatz (|C1|) transforms to be 



(#14) Z 2 [0,n] 1 : g t 

Vcr = T)crT = -1 
V<rC e = VC 6 T = -r)C e 



IT 3 , 9C 6 = T°, 



(C9) 



.9<r = .9c 6 = it 3 



(s,j/,«|x,j/,u) = -(x,y,u[x,y,w) = i(-l) x xT°,(C4) 
(x + l,y,w\x,y,u) = (x,y + l,v\x,y,u) =-i X T°, 
-(x,y,v[x,y,w) = (x + 1, y - 1, w[x, y, v) = xt 3 . 

Of course they belong to the 20 Z 2 spin liquids summa- 
The projected symmetry group (PSG) corresponds to rized in TABLE [TT1 



(#16) Z 2 [0,tt]S 

V<r = Va-T = -1, 
9<tC 6 = VC 6 T = -VC e = 1 



(CIO) 
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3. Four possible Zi spin liquids around (7(1) 
SL-[0, 7r] state: mean-field ansatz 

a. Consistent conditions on mean-field bonds 

Implementing the generic conditions mentioned ear- 
lier on several near neighbor mean-field bonds with PSG 
(|C6j) - (|C10[) . we obtain the following consistent condi- 
tions: 

(0) For on-site chemical potential terms A. s (x, y, s) — 
\{x,y,s) ■ r, translations operations Gt x 3 2i,2 m PSG 
guarantee that A s (x, y, s) = A s (0, 0, s) = A s , s = u,v, w. 
They satisfy 

g T A s g f T = -A.; (Cll) 
9crA u gX = A„, g a A v g\, = A w , g<,A w g\ = A u ; 
9c 6 A u g f C6 = A v , {gc e T 3 )A v (g C(i T 3 y = A Wi 
go e A w g f Ce = A u . 

(1) For 1st neighbor mean-field bond u a = [0,0, v}^ 
(there is only one independent mean-field bond, meaning 
all other 1st neighbor bonds can be generated from [0,0, v] 
through symmetry operations) 

g TU y T = ul (C12) 
{g<r9 2 c T 3 )ul{ gcT g 2 c j^ = -u a . 

(II) For 2nd neighbor mean- field bond Ub = [0, 1, w] we 
have 

ffruftflr = Ub, (C13) 
(ga-gc e )ub(ga-gc e ) 1 = -u\. 

(II) For 3rd neighbor mean-field bonds u c i = [1,0, u] 
and u c2 = [1, — 1, u] we have 

gTU c ig f T = -u c x, (C14) 
(ffc^HiG?^ 3 ) 1 = "Ii- 
and 

gTU c 29 f T = -U c2 , (C15) 

9<rU c2 gl = u\ 2 , 

(5o 8 T 3 )«e2(ffc 8 T 3 ) t - -<4- 

6. Mean-field ansatz of the four Z2 spin liquids near U(l) 
SL-[0,w] state 

For Z 2 [0,ir]a state with o CT = gc 6 = t° the mean-field 
ansatz are (up to 3rd neighbor mean-field bonds) 

u a = iaoT + o-it 1 , Ub — iboT°, (C16) 

Ucl = C3T 3 , U c2 = C 2 T 2 , 

A a = A3T 3 , s = u,v, w. 



Since we are considering a phase perturbed from the 
U(l) SL-[0,7r] state, we shall always assume a 7^ 
(1st neighbor hopping terms) in the following discussion. 
A Z 2 [0,ir]a spin liquid can be realized by 1st neighbor 
mean-field singlet pairing terms with a\ 7^ 0. 

For Z 2 [0, 7r]/3 state with g a = r°, gc e = it 3 the mean- 
field ansatz are (up to 3rd neighbor mean-field bonds) 

u a = i»oT° + air 1 , u b = ib T° + b\T X , (C17) 

Ucl = c 2 t 2 + c 3 r 3 , u c2 = 0, 
A„ = A 2 r 2 + A 3 r 3 , A v>w = -A 2 t 2 + A 3 t 3 . 

A Z 2 [0,tt]/3 spin liquid can be realized by 2nd neighbor 
pairing terms with a$bi — aibg =/= 0. 

For i?2[0,7r]7 state with g CT = ir 3 , gc 6 = t° the mean- 
field ansatz are (up to 3rd neighbor mean-field bonds) 

u a = ia r°, u b = i6 r° + biT 1 , (C18) 
Ucl = c 3 r 3 , u c2 = 0, 
A s = A 3 r 3 , s = u,v,w. 

A Z 2 [0,n]"/ spin liquid can be realized by 2nd neighbor 
pairing terms with b± 7^ 0. 

For Z 2 [Q, tt]S state with g a = gc 6 = It 3 the mean-field 
ansatz are (up to 3rd neighbor mean-field bonds) 

u a = ia T°, u b = ib T°, (C19) 
Uci = c 2 t 2 + c 3 r 3 , u c2 = 0, 
A s = A 3 r 3 , s = u,v,w. 

A Z2[0,7r](S spin liquid can be realized by 3rd neighbor 
pairing terms with c 2 =/= 0. 

4. Low-energy effective theory 

The reciprocal unit vectors (corresponding to unit vec- 
tors 04,2) on a kagome lattice are 61 = ^(x — -^y) and 
b 2 = i-^y, satisfying di ■ bj ~ Sij. In the mean-field 
ansatz (|U4)l of U(l) SL-[0,7r] the unit cell is doubled 
whose translation unit vectors are A\ = 2a\ and A 2 = a 2 . 
Accordingly the 1st BZ for such a mean- field ansatz is 
only half of the original 1st BZ with new reciprocal unit 
vectors being B\ = bi/2 and B 2 = b 2 . Denoting the mo- 
mentum as k = {k x ,k y )/a — kiBi + k 2 B 2 with |fci,2| < tt, 
we have 

fci = 2k x , k 2 = (k x + Vsk y )/2. (C20) 

The two Dirac cones in the spectra of U(l) SL-[0, it] state 
(|C4|) are located at ±Q with 

Q = (0,-|) = p 2 (C21) 

with the proper chemical potential A(i) = — x(V3 — 
l)r 3 added to mean-field ansatz (|C4[) . 
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For convenience we choose the following basis for 
Dirac-like Hamiltonian obtained from expansion around 
±Q: 

<t>+,t,A = Te^^ ■ 
(e-' 1 ^, 0, e'ra^, 0, 0, 0, e" 1 A^, 0, e' 1 A"", 0, y/2, 0) T , 

</>+,t,b = 7f e_iil7r • 
(1,0, e-' 1 !*, 0, v^e" 1 A", 0,-1, 0, e" 1 ^, 0, 0, 0) T , 
<P-.t,b = RrAfa = 0, k 2 = -f )<£+ ,t,6, 

</>±4,b = R r4>±,\.b- (C22) 

where ± are valley index for two Dirac cones at ±Q 
with Pauli matrices /x and b = A, B are band indices 
(for the two bands forming the Dirac cone) with Pauli 
matrices v. Pseudospin indices £ =f , 1 are assigned to 
the two degenerate bands related by time reversal, with 
Pauli matrices cr. The corresponding creation operators 
for these modes are "f^ s b = "0±q0±,s.6 in the order 
of (0, 0, u), (0, 0, v), (0, 0, to), (1, 0, u), (1, 0, u), (1, 0, w) for 
the six sites per doubled new unit cell. Notice that in 
terms of /-spinons we have ip^ = (fl, fi). 

-1 0" 

10® gr, i?T 2 (k) = 
1 



Here Rt 



'2x2 



'6x6 



g T2 and R Tl (k) 



- i ki 



o 



t/Ti are transformation matrices on 12- 



1 
0-10 
0-1 

component eigenvectors for time reversal T and transla- 
tion T\,2 operations. By definition of PSG the eigenvec- 
tors 0k with momentum k = k\B\ + k 2 B 2 = (fci, k 2 ) and 
energy E have the following symmetric properties: 



(kiM) 



T 2 : 4>(kiM)( E ) = R TA k iM)(t>(kiM)( E )- 

cj> and 4> are the basis after and before the symmetry 
operations. 

In such a set of basis the Dirac Hamiltonian obtained 
by expanding the U(l) SL-[0, tt] mean-field ansatz (|C4[) 
around the two cones at ±Q is 

H miac = -^^VC-fei/ 1 + fe^ 2 )* k (C23) 
k v 2 

k should be understood as small momenta measured 
from ±Q. Possible mass terms are 

^0,1,2,3^1,2^0 

and 



/' 



0,1,2,3 0,3^3 



However not all of them are allowed by 



symmetry. Here we numerate all symmetry operations 
and associated operator transformations: 
Spin rotation along z-axis by angle 6: 

i 4. ; e 



Spin rotation along j/-axis by tt: 

*k -> * T kM 2 c r2 ^ 2 

Time reversal T: 

* k ^*k(-^ 2 ) 

Translation T\: 
Translation T 2 : 

Considering the above conditions, the only symmetry- 
allowed mass terms are ^^mi^k with mi = //V 1 ^ 
and 777,2 = p?o i v i . 

The transformation rules for mirror reflection er and 
7r/3 rotation C6 depend on the choice of g^gce m the 
PSG. In general we have 

a : ^ -> < k A/ CT ( 5(T ), 
Ce: * k ^*c 6 k M c 6 (ffc 6 ). 



Using the basis (|C22[) the 8x8 matrices M^Cg can be 
expressed in terms of Pauli matrices /x <g) er <g) f . For the 
four spin liquid we have 



M CT ( 5(T = r°) = /j 3 ® cr° (g 
M a (g a = ir 3 ) = ^®a 3 

M C8 (ffCe =T°) 

Mc 6 (5c 6 - ir 3 ) 



1 
i 



i 

-1 



e-'ra' 
; -iA w 

e i A 7r 
e'A^ 



It turns out in Z 2 [0,ir](3 state, only the 1st mass term 
?77. 1 = /i'W is invariant under er and C§ operations. 
In other 3 states neither mass terms 7711,2 are symmetry- 
allowed. As a result we only have one gapped Z 2 spin 
liquid, i.e. Z 2 [0,ir]ft state in the neighborhood of U(l) 
Dirac SL-[0,7r] state. 

Let's consider mean-field bonds up to 2nd neighbor for 
ansatz Z 2 [0,ir)(3. Perturbations to the two Dirac cones 
of U(l) SL-[0, 7r] with A3 = (v3 — l)cso in general has the 
following form 

SH = [A 3 - (\/3 - l)oo - (\/3 + l)6o] yPo 3 ^ 
+ [(^3 + l)6i - A a - (a/3 - l)ai]/iVi/° (C24) 

This means we need either 1st neighbor (<xi) or 2nd neigh- 
bor (&i) pairing term to open up a gap in the spectrum. 
Meanwhile these pairing terms break the original U{1) 
symmetry down to Z 2 symmetry. 
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Appendix D: Z-2 spin liquids in the neighborhood of 
uniform RVB state 

The mean-field ansatz of the uniform RVB state is sim- 
ple: 

H M F = X E ft<rh° ( D1 ) 

<ij>,<7 

where x is a rea l parameter and < ij > represents 
sites i,j being nearest neighbor (n.n.) of each other. 
It's straightforward to show the PSG elements of such 
a mean-field ansatz are 

G u {x,y,s)=g u , U = T 1>2 , T, a, C 6 . (D2) 

and SU(2) matrices gu satisfy 

ffT r 3 . g ; = -r 3 , (D3) 
9uT 3 glj = T 3 , U = T h2 ,<J,C 6 . 

It turns out there are only 4 gauge-inequivalent Z 2 PSGs 
as solutions to (Ta1)) - (|A16|) with the form (fD2|) . In other 
words, there are only 4 different Z 2 in the neighborhood 



of a uniform RVB states. Choosing a proper gauge they 
all satisfy gx = It 1 , gr i2 = T ° an d Vt 12 t = ?7i2 = 
VCgT! 2 = Vo-t 12 = lj Vt = —1- These four states are 
characterized by: 

(#1) Z 2 [Q,0}A: grT =g Ce =r a , (D4) 

VctT = VC e T = V<r = VC 6 = VcrC 6 = 1- 

(#5) Z 2 [0,0}B : g a = t°, 5C 6 = ir 3 , (D5) 

= 7?o- = 1, 7?C 6 X = ??C 6 = '?<tC„ = -1. 

(#15) Z 2 [0,0]C : 5ct = ir 3 , g Ce = r°, (D6) 

= 7?o- = -1, 7/C 6 T = ??C 6 = f? CT C 6 = 1- 

(#13) Z 2 [0,0}D : , 9ct = g Ce = ir 3 , (D7) 

»? CT T = VCeT = Vo- = VC 6 = V<rC 6 = ~1- 

It turns out these four SLs around uniform RVB 
state are all gapped as shown in TABLE [TTJ 
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